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Abstract 



We present a new view of Feynman diagrams for the field theory of 
transport on percolation clusters. The diagrams for random resistor net- 
works are interpreted as being resistor networks themselves. This simpli- 
fies the field theory considerably as we demonstrate by calculating the frac- 
tal dimension Db of the percolation backbone to three loop order. Using 
renormalization group methods we obtain D B = 2 + e/21 - 172e 2 /9261 + 
2e 3 (-74639 + 22680C (3)) /4084101, where e = 6 - d with d being the spatial 
dimension and C (3) = 1.202057... 

PACS numbers: 64.60.Ak, 64.60.Fr, 72.80.Ng, 05.70.Jk 
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Percolation has gained a vast amount of interest over the last decades (for a review see 
e.g. UIU). Though it represents the simplest model of a disordered system it has many 
applications e.g. polymerization, porous and amorphous materials, thin films, spreading of 
epidemics etc. Consider a <i-dimensional lattice where each bond is randomly occupied 
with probability p or empty with probability 1 — p. Occupied and empty bonds may stand 
for different physical properties. Assume that occupied bonds are electrical conductors 
whereas empty sites are insulators and that currents can flow only between nearest neighbors. 
Suppose a potential difference is applied between two sites x and x' located on the same 
cluster. In general not all bonds do carry non-zero current since there may be dangling 
ends. This gives rise to the notion of the backbone. It is defined as the set of bonds 
which are connected to both x and x' by mutually non-intersecting paths. Except for 
Wheatstone bridge type configurations these are the bonds that carry non-zero current. 
The fractal dimension Db of the backbone is defined near the critical concentration p c by 
Mb ~ \x—x'\ Db , where Mb denotes the average number of bonds (the mass) of the backbone. 

In this paper we evaluate Db by renormalized field theory. Our approach is based on a 
field theoretic formulation of the randomly diluted nonlinear resistor network by Harris || 
which itself was based on work by Stephen [|J and Harris and Lubensky ||. The aim of this 
letter is to present our new interpretation of Feynman diagrams as being resistor networks 
themselves and to employ this interpretation to derive Db up to third order in e = 6 — d. 

Consider a nonlinear generalization of the random resistor network as proposed by Kenkel 
and Straley . The bonds between sites i and j obey a generalized Ohm's Law 

\'j V t /',,/„/,,/ ' (1) 

or equivalently 

0ij{Vj-V i )\V s -V i \- x = I i j 1 (2) 

where cjjj (pij) is the nonlinear conductance (resistance) of the bond, ij-j is the current 
flowing through the bond and Vi is the potential at site i. The exponents r and s are 
describing the nonlinearity with r = s^ 1 . 

The power P dissipated on the backbone between x and x' of this nonlinear network 
reads 

P = UI = R r (x,x')\I\ r+1 , (3) 

where U denotes the voltage between the two ports, I the resulting current and R r (x,x') 
the resistance of the backbone. On the other hand we may write 

p = T,\ v j- v i\hj = T,PiMA r+1 , (4) 

where the sum is taken over all bonds on the cluster. The limit r — > — 1, taken from above, 
provides for a convenient way of summing up all conductors carrying non-zero current: 

R-i(x,x') = ^Pi,j ■ (5) 

hi 
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We restrict ourself to the case that all conductors have identical resistance p. Hence is 
proportional to Mb and Db is identical to lim^.! <\> r j v, where v is the correlation length ex- 
ponent and r is the resistance exponent defined by M r = (x(x, x')R r (x, x'))c/ (x( x i x '))c ~ 
\x — x'\^ r ' v '. (...)c denotes the average over all configurations of the diluted lattice and 
x(x, x') is an indicator function that takes the value one if x and x' are on the same cluster 
and zero otherwise. 

The resistance R r (x,x') can be obtained by solving the circuit equations 

J2^Av l -v J )\v l -v,r 1 = i l , (6) 

j 

where Jj = I {8^ x — 5i jX >). The circuit equations may be viewed as a consequence of the 
variation principle 



_d_ 

dVi 



(7) 



where {V} denotes the set of voltages belonging to the sites of the backbone. Obviously 
the backbone may contain closed loops as sub-networks. Suppose there are currents {-f^j 
circulating independently around these closed loops. Then the power is not only a function 
of / but also of the set of loop currents. Conservation of charge holds for every ramification 
of the backbone and this gives rise to another variation principle: 

9 P({J(0},J)=0. (8) 



did) 



Eq. (H) may be used to eliminate the loop currents and thus provides us with a method to 
determine the total resistance of the backbone via Eq. (H). 

A field theory for the nonlinear random resistor network was set up by Harris || in 
analogy to the linear model |H|. In order to overcome difficulties associated with {...)c 
one employs the replica technique |J. The network is replicated .D-fold: V x — ► V x — 
(V^, . . . , V^ ^. One considers the correlation function G (x, x'; Xj = (tp^(x)tjj_^(x')^ of 

ip x (x) = exp (iX ■ V x ) where \-V x = Y, a \ {a) V^ and A ^ 0: 



g (x, x'- x) = [z-° j n n dv ? ex p {-^i p ({?}) + zX ■ - &)) ) • ( 9 ) 

Here P ({^}) = J2i,j,a ~~ Vj a \ s+1 and Z is the usual normalization. In contrast 

to the linear network P is not quadratic and hence the integration is not gaussian. As a 
working hypothesis we assume that a saddle point approximation is justified. For details 
and conditions to be imposed on A see ||. The saddle point equation is nothing else than 
the variation principle stated in Eq. (0). Thus the maximum of the integrand is determined 
by the solution of the circuit equations (§) and, up to an unimportant constant, 

G (x, x'- X) = ^exp (J^Rr & a0) ) = (*(*, *))c ( X + ^ M -( x ' x + ■ ■ ■) > ( 10 ) 
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n / i \2\(H"l)/2 

where A r = X)a=i ( ~~ ^ ) ■ Note that the limit D — > has to be taken before r — > — 1 
for Eq. fllCf ) to be well defined. Contact to the Potts-model can be established by switching 
to voltage variables 9 = A6k taking discrete values on a D-dimensional torus, i.e. k is chosen 
to be an /^-dimensional integer with — M < k^ < M and k^ = k^mod{2M). In this 
discrete picture there are (2M) D — 1 independent state variables per lattice site and one 
introduces the Potts-spins 

9 ff (x) = (2M)- D E exp (iX ■ 6) i^x) = 8^ - (2M)~ D (11) 

subject to the condition J2g&g( x ) — 0. 

The replication procedure induces the effective Hamiltonian 

^e P = -ln(exp(--l T p)) c (12) 

which may be expanded in terms of i/j: 

Hre P = - E E K ( X ) i>X (*) i>-X V) ■ (13) 
x > x ' A^O 

Next the kernel is Taylor expanded in the limit of large a: 

k(X) = t- wA r , (14) 



with t and w ~ a 1 being expansion coefficients and higher order terms are neglected since 
H rep is decaying exponentially. By defining the discrete derivative d/dO^ through 

-E**(*) = E A (a) ^*)tM*) (is) 

" A>0 

one obtains upon Fourier transformation 

tf(A 9 -)=r- m (A 9 -) (r+1,/ \ (16) 

To set up a field theoretic Hamiltonian W we proceed with the usual coarse graining step and 
replace the Potts-spins $^ (x) by the order parameter ip (x, 9^ defined on a <i-dimensional 
spatial continuum. Constructing all possible invariants of the symmetries of the model 
from J2q ¥ ( x ) 0) (p denotes some power > 1) and gradients thereof leads to the following 
Hamiltonian in spirit of the Landau-Ginzburg- Wilson functional (for details see ||): 

H = / AE{!/ + \ (v„)- - > (4rr*%+ V} , (17) 

where terms of higher order in the fields have been neglected since they turn out to be 
irrelevant in the renormalization group sense. Note that Tt reduces to the usual Potts- 
model Hamiltonian by setting w = 0. 
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Now we set up a diagrammatic expansion. Contributing elements are the vertex —g and 

> nrnnfl cmtnr 



the propagator 

l-^o 1 *X0 (lg) 



p 2 + r — toA r p 2 + r — wA r p 2 + r 

Eq. (]T8|) shows that the principal propagator decomposes into a propagator carrying A's 
(conducting) and one not carrying A's (insulating). This allows for a schematic decompo- 
sition of principal diagrams into sums of diagrams consisting of conducting and insulating 
propagators. Here a new interpretation of the Feynman diagrams emerges ||. They may be 
viewed as resistor networks themselves with conducting propagators corresponding to con- 
ductors and insulating propagators to open bonds. Schwinger parameters Sj of conducting 
propagators correspond to resistances a" 1 and the replica variables i\ to currents. The 
replica currents are conserved in each vertex and we may write A, = A, {X, {ft}), where A is 

an external current and {k} denotes the set of independent loop currents. The A-dependent 
part of a diagram can be expressed in terms of its power P: 



exp wVsjAri = exp(u>P(A, {£})) • (19) 




The new interpretation suggests an alternative way of computing the Feynman diagrams. 
To evaluate sums over independent loop currents 



J2exp(wP{X,{K})) (20) 
{«} 

we employ the saddle point method. Note that the saddle point equation is nothing else 
than the variation principle stated in Eq. Thus solving the saddle point equations is 
equivalent to determining the total resistance R({si}) of a diagram and the saddle point 
evaluation of (|20|) yields 

exp (Rr ({s,}) wA r ) . (21) 

A completion of squares in the momenta renders the momentum integrations straightfor- 
ward. Thereafter all diagrams are of the form 

/ (p 2 , A 2 ) = I P (p 2 ) + I w (p 2 ) wA r + . . . 

= f o °°T[ds i {l + R r ({s l })wA r + ...]D(p 2 ,{s l }) , (22) 

where D (p 2 , {s^}) is a usual integrand of the 3 -theory. The 3 -theory was investigated 
to three loop order by de Alcantara Bonfim et al |T3] and hence the remaining task is to 
calculate the contributions proportional to w. 

In order to check if our working hypothesis holds we performed two loop calculations for 
the cases r — > and r — > oo and compared to known results. In the limit r — > the resistance 
between two points becomes essentially equal to the length of the shortest paths between 
these points. We mapped our diagrams onto those studied by Janssen and obtained ex- 
actly the same diagrammatic expansion. Consequently, our result for the exponent governing 
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the so-called chemical distance d min = 2 - e/6 - [937/588 + 45/49 (In 2 - 9/10 In 3)] (e/6) 2 + 
O (e 3 ), e = 6 — d, is the same as given in ||. The limit r — > oo is related to the red (singly 
connected) bonds. Our calculation gives unity for the corresponding exponent in accordance 
with results by Blumenfeld and Aharony |HJ and de Arcangelis et al. [11 1. We rate these 
two loop results as a strong indication for the validity of the saddle point approach. 

Now we turn to the calculation of Db- In the limit r —>■ —1 only the non-planar diagrams 
listed in Fig. [I] contribute to the diagrammatic expansion. We use dimensional regularization 
and renormalize w — > Z~ x Z w w. By employing minimal subtraction to compensate e-poles 
we obtain 



u 2 u 3 
4e e 2 



!-*U-H«s) e 

12 144 3 . 



+ o 



(23) 



where u oc g 2 fi~ e with \i being an inverse length scale. Note that all non-primitive divergen- 
cies are cancelled as renormalizability of the perturbation expansion requires. The critical 
exponents are determined by the Wilson-functions 7... (u) = fi-^-lnZ,., evaluated at the in- 
frared stable fixed point u* . In particular we are interested in rj = 7 (u*) and ip = j w (u*) 
governing the scaling relation 



G (|x - x'|; w) = l d ~ 2+n G (Z|x - x'|; w/l 2 ~^) 



where / is a inverse length scale, r\ was calculated to order e in [12|. For ip we find 



if, 



+ 



16C(3) 



2075 



From Eq. 



and hence 



in conjunction with Eq. (|T 



126 

it follows that 



C(c*) 



D 



B 



G 



2-r] + tp 



w) 



1 

2 + — e 
21 



_/|2— tj-H0 



+ 



172 



-74639 + 22680C (3) 



(24) 



(25) 



(26) 



(27) 



9261 4084101 
Note that our result agrees to second order in e with calculations by Harris and Lubensky 
fH| based on another approach. This is again in favor of our working hypothesis. 

We compare our result to numerical simulations by Grassberger |14| and Moukarzel [15 



Due to the rich structure of 77 in the percolation problem ip is better suited for such a 
comparison than Db- It is known exactly that ip vanishes in one dimension. This feature is 
incorporated by rational approximation yielding 

2e 2 



V « -— ( 1 - - ) ( 1 + 1.2625 



500 



(28) 



which is compared to simulations in Fig. 0. For d = 4 the results agree within the numerical 
errors. However, a higher accuracy of the numerical estimate is desirable. For d = 3 and 
d = 2 the analytic result looks less realistic and the numerical values are larger. The shape 
of the dependence of ip on dimensionality is much the same. 

We conclude with a few comments. Our interpretation of the Feynman diagrams simpli- 
fies calculations considerably. The technique used here can be applied to study other aspects 
of transport on percolating clusters. In d — 4 our result for Db agrees with recent numerical 
simulations. For dimensions close to the upper critical dimension six, our result is the most 
accurate analytical estimate for Db that we know of. 
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FIGURES 



FIG. | 

The diagrams we computed to determine D B . 
the solid dots for l^-insertions. 



The lines stand for conducting propagators, 



FIG. | 

Dependence of the exponent ip on dimensionality. The rational approximation (triangles) is 
compared to numerical results (circles) by Grassberger (d = 2) and Moukarzel (d = 3,4). 
They determined Db = 2 — r] + ip = r )/v + il)hY simulations. For d = 2 we insert the exact 
TT6t|T7f v = 4/3 and 7 = 43/18. For d = 3 we use Monte Carlo results by Ziff and 
= 0.875 ± 0.008, 7 = 1.795 ± 0.005. For d = 4 we take v' 1 



values 
Stell 



and 7 



: v - 
1.44 



43/18. For d 
= 1.795 ± 0.005. For d 



1.44 ±0.05 15 
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FIG. 1. 
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e = 6-d 



FIG. 2. 
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